Draft version December 6, 2011 

Preprint typeset using DTjrpC style emulateapj v. 11/26/03 


STELLAR POPULATIONS OF HIGHLY MAGNIFIED LENSED GALAXIES: YOUNG STARBURSTS AT Z ~ 2. 

Eva Wuyts 1 ' 2 , Jane R. Rigby 3 , Michael D. Gladders 1 ’ 2 , David G. Gilbank 4 , Keren Sharon 2 , Megan B. Gralla 1,2 , 

Matthew B. Bayliss 1,2 

Draft version December 6, 2011 



CNf 


o 

QJ 


IT) 




T 

Q 


oo 

m 


CT 

m 

cn 

oo 

CN 


N*****’ 





ABSTRACT 

We present a comprehensive analysis of the rest- frame UV to near-IR spectral energy distributions 
and rest-frame optical spectra of four of the brightest gravitationally lensed galaxies in the literature: 
RCSGA 032727-132609 at 2 - 1.70, MS1512-cB58 at 2 - 2.73, SGAS J152745. 1+0652 19 at * = 2.76 
and SGAS J122651. 3+215220 at 2 = 2.92. This includes new Spitzer imaging for RCSGA0327 as well 
as new spectra, nearTR imaging and Spitzer imaging for SGAS1527 and SGAS1226. Leasing magni- 
fications of 3-4 magnitudes allow a detailed study of the stellar populations and physical conditions. 
We compare star formation rates as measured from the SED fit, the Iia and [O II] A3 72 7 emission 
lines, and the UV+IR bolometric luminosity where 24 pm photometry is available. The SFR estimate 
from the SED fit is consistently higher than the other indicators, which suggests that the Calzetti dust 
extinction law used in the SED fitting is too flat for young star- forming galaxies at z ~ 2. Our analysis 
finds similar stellar population parameters for all four lensed galaxies: stellar masses 3 — 7 x 1O 9 M 0 , 
young ages ~ 100 Myr, little dust content E(B — V)=0. 10-0.25, and star formation rates around 20- 
100 M 0 yr“U Compared to typical values for the galaxy population at z ~ 2, this suggests we are 
looking at newly formed, starbursting systems that have only recently started the build-up of stellar 
mass. These results constitute the first detailed, uniform analysis of a sample of the growing number 
of strongly lensed galaxies known at z ~ 2. 

Subject headings: galaxies: high-redshift, strong gravitational lensing, infrared: galaxies 


1. INTRODUCTION 

Recent years have seen a significant rise in the number 
of known gravitationally lensed galaxies at z = 1 — 3, 
both reported from well-defined survey searches (Bolton 
et ah 2006, Cabanac et al. 2007, Hennawi et al. 2008, 
Lin et al. 2009, Gladders et al. 2011 in preparation) 
and from serendipitous discoveries (Allam et al. 2007; 
Belokurov et al. 2007; Smail et al. 2007). The lens- 
ing magnification grants access to high signal- to- noise 
photometric and spectroscopic observations over a wide 
range of wavelengths and creates unique opportunities 
for very detailed study of the stellar populations, dy- 
namics and physical conditions of star- forming galaxies 
at 0 ~ 2, complementing statistical studies of large sam- 
ples of unlensed star-forming galaxies at these redshifts 
(e.g. Shapley et al. 2005, Erb et al. 2006a, Daddi et 
al. 2007, Reddy et al. 2010). Today, these studies are 
limited by systematic uncertainties related to the recipes 
used to estimate physical parameters such as metallic- 
ity, dust extinction, star formation rate and star forma- 
tion history. Until the era of extremely large telescopes, 
lensed galaxies provide unique opportunities to increase 
our physical understanding of the composition and as- 
sembly history of stellar populations during this period 
of peak star formation in the Universe. 

* Based in part on observations collected at the 3.5 m Apache 
Point Observatory telescope in New Mexico, which is owned and 
operated by the Astrophysical Research Consortium 
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So far most lensing systems have been studied indi- 
vidually, based on different follow-up observations and 
analysis techniques. To take advantage of the growing 
number of known strongly lensed sources, it is crucial to 
invest in a uniform study of their properties. In this pa- 
per we present a full analysis of multi-wavelength pho- 
tometry and rest-frame optical spectroscopy of four of 
the brightest distant lensed galaxies known to date. RC- 
SGA 032727-132609 was recently discovered in the Sec- 
ond Red-Sequence Cluster Survey (Gilbank et al. 2011); 
initial analysis of optical/near-IR photometry and con- 
struction of a lens model for the foreground galaxy clus- 
ter are described in Wuyts et al. (2010). Rest-frame op- 
tical spectroscopy from Keek/NIRSPEC is presented in 
Rigby et al. (2011). This paper adds Spitzer imaging 
to the SED to probe the rest- frame near-IR. MS 1512- 
cB58 has been studied extensively since its discovery 
(Yee et al. 1996) and a multitude of photometric and 
spectroscopic data is available for this system (Elling- 
son et al. 1996: Teplitz et al. 2000; Siana et al. 2008). 
SGAS J 152745. 1+0652 19 and SGAS J 12265 1.3+2 15220 
are part of the SDSS Giant Arcs Survey (Bayliss et al. 
2011); optical photometry and lens models are reported 
in Koester et al. (2010), Here we present near-IR and 
Spitzer imaging as well as rest-frame optical spectra for 
both sources. We will refer to these galaxies as RC- 
SGA0327, cB58, SGAS1527 and SGAS1226. 

The paper is organized as follows. §2 describes the 
near-IR and Spitzer imaging and the rest- frame optical 
spectroscopy. The methods used for the photometric 
analysis and spectral energy distribution modeling, as 
well as the various reddening and star formation rate di- 
agnostics are discussed in §3. §4 presents the physical 
conditions and stellar populations of the lensed galaxies 
as derived from these methods. We adopt a flat cosmol- 
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ogv with Qm = 0.3 and H 0 — 70kms 1 Mpc x . Ail 
magnitudes are quoted in the AB system. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Imaging 

Optical and near-IR photometry of RCSGA0327 and 
optical photometry of SGAS1527 and SGAS1226 are 
taken from their respective discovery papers (Wuyts et 
ah 2010; Koester et ah 2010). For cB58, we take opti- 
cal and near-IR photometry from Ellingson et ah (1996) 
and Spitzer photometry from Siana et al. (2008). Here 
we present new Spitzer imaging for RCSGA0327 as well 
as new T near-IR and Spitzer imaging for SGAS1527 and 
SGAS1226. 

We observed SGAS1527 and SGAS1226 with the 3.5 m 
telescope at the Apache Point Observatory (APO) in 
New Mexico in the J-, H- and K s - bands during four 
half nights on 2010 January 30, February 24, March 4 
and March 29 with seeing of 077-079. A pipeline of stan- 
dard IRAF 1 tasks is used to sky-subtract and stack the 
dithered images, giving total integration times of 3600 s 

- 4500 s - 3600 s in J - II - K s for SGAS1527 and 3600 s 

- 3150 s - 2700 s in J - H - K s for SGAS1226. 

IRAC (Fazio et al. 2004) and MIPS (Rieke et ah 2004) 
observations of RCSGA0327 and SGAS1527 were ob- 
tained through Spitzer program 50823 (PI: M. D. Glad- 
dens). RCSGA0327 was observed for 600 s in all four 
IRAC bands (3.6/4.5/5.8/8.0/im) and 700s in MIPS 
24 /xm; SGAS1527 was observed for 3000 s in all IRAC 
bands and 1300 s at 24 pm. 600 s warm IRAC obser- 
vations of SGAS1226 at 3.6 and 4.5 pm were obtained 
through Spitzer program 70154 (PI: M. D. Gladders). 
The data is reduced with the MOPEX software dis- 
tributed by the Spitzer Science Center (SCC) and driz- 
zled to a finer pixel scale of 075 pix” 1 for IRAC and 
175 pix” 1 for MIPS. 

2.2. Rest- frame optical spectroscopy 

Keck/NIRSPEC rest-frame optical spectra were pub- 
lished for cB58 by Teplitz et al. (2000) and for RC- 
SGA0327 by Rigby et al. (2011). Here we present new 
spectra for SG AS 1527 and S GAS 1226. 

We observed SGAS1527 and SGAS1226 with the NIR- 
SPEC spectrograph (McLean et ai. 1998) on the Keck II 
telescope on the night of 2010, February 4. The weather 
was clear, and the seeing was measured as 0.85" and 
0.45" when the telescope was focused during the night. 
We used the low- resolution mode and the 0.76" x 42" 
slit which was rotated to position angles of 112° and 
121° East of North, respectively, to maximize coverage 
of each arc. Finder charts showing the location of the slit 
on each target are given in Figure 1. The targets were ac- 
quired by offsetting from the brightest cluster galaxy on 
the near-IR slit-viewing camera; target acquisition was 
verified by direct imaging with this camera. The targets 
were nodded along the slit in an AB pattern, with expo- 
sures of 600 s. Table 1 summarizes the filters and total 

1 IRAF (Image Reduction and Analysis Facility) is distributed 

by the National Optical Astronomy Observatories, which are op- 
erated by AURA, Inc., under cooperative agreement with the Na- 
tional Science Foundation. 


Table l. Spectroscopic Observation Log. 


source 

filter 

time 

s 

wavelength range 
fim 

SGAS1527 

NIRS PEC-4 

1800 

1.302-1.585 


NIRSPEC-6 

3000 

1.770—2.190 


NIRS PEC- 7 

1620 

2.242-2.521 

SGAS1226 

NIRS PEC- 5 b 

1800 

1.447-1.730 


NIRS PEC- 7b 

3000 

1,864-2.285 


integration times. For SGAS1527, the default grating 
angles were used, but for S GAS 1226 the grating angles 
were adjusted to give the wavelength ranges reported. 
The A0V stars 7 Ser and HD 109055 were observed every 
hour as telluric standards for SGAS1527 and SGAS1226 
respectively. In addition, V. Tilvi and J. Rhoads of the 
Arizona State University have kindly obtained additional 
spectra of SG AS 1527 on 2010, September 17, with the 
same setup. The AOV star HD 12021, observed 10 hr 
later, is used as a telluric standard. Due to the time de- 
lay, the absolute fluxing in this bandpass should not be 
trusted. 

The spectra are reduced with the nirspemreduce pack- 
age written by G. D. Becker. The data reduction, flux- 
ing, and line-fitting procedures are described in Rigby et 
al. (2011). Solar abundances are taken from Table 1 of 
Asplund et al (2009). The reduced spectra are plotted 
in Figures 2 and 3; line fluxes are reported in Table 2. 
At observed wavelengths where the Earth's atmosphere 
has very little transmission, small mismatches in airmass 
and precipitable water vapor (PWV) between the tel- 
luric and the science target can have a large effect on the 
flux calibration. To quantify this, Table 2 reports the ex- 
pected atmospheric transmission 2 for each detected emis- 
sion line integrated over the gaussian line width as given 
in §4.1. The quoted uncertainty combines the measure- 
ment uncertainty in the line widths, the systematic un- 
certainty in PWV and the systematic uncertainty from 
a 0.1 airmass mismatch between the telluric and science 
target. Subsequent calculations involving line fluxes will 
take into account both the la flux uncertainties and the 
fractional uncertainty on the atmospheric transmission. 

3. METHODS 
3.1. Photometry 

The photometric analysis of the imaging data follows 
the method detailed in Wuyts et al. (2010). In short, the 
data are transformed to a common reference and pixel 
scale and empirical, normalized point spread functions 
(PSF) are created for each image. We define object aper- 
tures by tracing a curve along the extended source and 
convolving it with the appropriate PSF. A series of aper- 
tures of increasing radial extent are defined as isophotes 
of this convolution. The apertures can be very much 
non- circular and are described by equivalent radii based 
on circular apertures that extend to the same isophotes. 

2 The atmospheric transmission is baaed on the model absorp- 
tion spectra for Mauna Kea published by the Gemini Observatory: 
http: //www .gemini.edu /sciops/ telescopes-and-sites/ observing- 
condlt ion-const raints/ir- transmission-spectra. The PWV at the 
time of observations was derived from historical atmospheric 
opacity data at http://puuoo.caltech.edu/. 



Fig. 1. Finder charts showing NIRSPEC slit positions for SGAS1527 (left) and SGAS1226 (right), both 45 x 45". The CFHT Megacam 
r-band image of SGAS1527 and the Gemini CMOS t-band of SGAS1226 (Koester et al. 2010) are overplotted with the NIRSPEC 0.76" x 
42" longslit at position angles of 112° and 121° East of North respectively. The missing region in the i-band image of SGAS1226 is due to 
the GMOS chip gap. 
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Fig. 2. NIRSPEC spectra for SGAS1527. The fluxed spectra are plotted in black, with the la error spectrum in red. The X-axis shows 
observed wavelength in Angstroms; the Y-axis shows observed flux in units of erg s""" 1 cm -2 . Blue labels mark detected emission lines as 
well as lines for which we measure upper limits. 


Table 2. Measured line fluxes. 


line ID 

^ obs 

/im 10"" 

flux 

16 erg s' 1 cm" 2 

atm trans 


SGAS1527 

- Filter N4 


[Oil] 3727 

1.40208 

114:1 

0.544:0.08 

[NelllJ 3869 

1.45578 

4±1 

0.624:0.07 


SGAS1527 

- Filter N6 


H 8 

1.82932 

9 ±2 

0.174:0.08 

[OIII] 4959 

1.86616 

20.0T0.2 

0.234:0.08 

[OIII] 5007 

1.88397 

18.04:0,3 

0.64T0.08 


SGAS1527 

- Filter N7 


Ha 

2.46925 

10±3 

0.824:0.03 


SGAS1226 ■ 

- Filter N5b 


[Oil] 3727 

— 

9.64:0.7 

0.76T0.06 


SGAS1226 • 

- Filter N7b 


m 

1.90895 

1.84:0.4 

0.554:0.09 

[OIII] 4959 

1.94697 

4.5±0.7 

0.744:0.04 

[OIII] 5007 

1.96641 

134:4 

0.674:0.03 


We use the GALFIT package (version 3.0, Peng et al. 
2010) to fit Sersic profiles to neighboring galaxies that fall 
within the object apertures and subtract these from the 
images. Accurate masking of neighboring galaxies be- 
comes very important for the Spitzer data due to the ex- 
tended, non-circular PSFs of the IRAC and MIPS instru- 


ments. We have found the most robust masking method 
to be based on cross-convolution. This consists of con- 
volving the Spitzer data with the PSF of a chosen optical 
band and similarly convolving the optical GALFIT mod- 
els of neighboring sources with the relevant Spitzer PSF. 
These convolved models are scaled as needed and sub- 
tracted from the convolved Spitzer image. 

We measure magnitudes at an equivalent radius of 
twice the FWHM of the image. The optical and near- 
IR data are aperture-corrected to an equivalent radius 
of 6" based on the curve of growth of the PSF reference 
star. In Wuyts et al. (2010), we convolved all optical 
and near~IR images of RCSGA0327 to the spatial res- 
olution of the //-band, which presents the worst seeing 
over the 9 optical/ near- IR bands, in order to measure 
flux from the same physical region of the source in all 
bands. This is not repeated here; by aperture correct- 
ing final magnitudes to an equivalent radius of 6", the 
total source flux in each band is measured, regardless of 
prior PSF- matching. We have reanalyzed the data for 
RCSGA0327 without matching the PSFs and find final 
magnitudes consistent with the photometry reported in 
Wuyts et al. (2010), within the la uncertainties. The 
aperture corrections for the IRAC data are applied as 
detailed in the IRAC handbook 3 , the correction for the 
MIPS data is calculated from a Tiny Tim model of the 
24 pm PSF. The MIPS photometry also requires a color 
correction of -0.04 mag, as stated in the MIPS handbook 4 
(assuming f v ~ v~ 2 ). We derive limiting magnitudes for 
a few non-detections in the Spitzer data: the counter- 
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Fig. 3. NIRSPEC spectra for SGAS1226. The fluxed spectra are plotted in black, with the 1 <j error spectrum in red. The X-axis shows 
observed wavelength in Angstroms; the Y-axis shows observed flux in units of erg s'" 1 cm” 2 . Blue labels mark detected emission lines as 
well as lines for which we measure upper limits. 


image of RCSGA0327 is not detected at 5.8 and 8.0 /i m, 
and SG AS 1527 is not detected at 24 //m. We convolve 
the image of the source in a chosen optical band with 
the relevant Spitzer PSF, scale it and manually insert 
it into the convolved Spitzer image at 100 random posi- 
tions away from real sources. Limiting magnitudes are 
determined from the scaling where SExtractor (Bert in & 
Arnouts 1996) detects 90% of these mock objects at 3 a. 

Final magnitudes are corrected for galactic extinction 
(Schlegel 1998). Photometric uncertainties include Pois- 
son noise, absolute zeropoint uncertainties, uncertainties 
from the Spitzer aperture corrections as detailed in the 
instrument handbooks, and uncertainties from the mask- 
ing of neighboring galaxies; the latter two dominate the 
total uncertainty. We independently analyze the IRAC 
and MIPS data for cB58 and find final magnitudes con- 
sistent with the photometry in Siana et al. (2008). We 
choose to use uncertainties twice as large as reported 
there, to take into account the influence of the halo of 
the neighboring cD galaxy. 

3.2. Spectral energy distribution modeling 

The available photometry fully characterizes the rest- 
frame UV to near-IR spectral energy distributions of 
these galaxies, allowing us to constrain the stellar popu- 
lations (stellar mass, age, extinction and star formation 
rate) through a comparison to stellar population synthe- 
sis models. The presence of strong emission lines at rest- 
frame optical wavelengths could affect the broad-band 
photometry and therefore the SED fit. From the emis- 
sion line fluxes derived from the rest-frame optical spec- 
troscopy, we estimate the contribution of the strongest 
emission lines to the H - and JFG-band magnitudes to be 
at most 5-10%, which does not have a significant effect 
compared to the estimated photometric errors. Reddy 
et ah (2010) have also shown that the effect of emis- 
sion lines on the best-fit age and stellar mass is further 
reduced when IRAC data (which is unaffected by line 
contamination at these redshifts) is incorporated. 

The fitting procedure, as well as its caveats and de- 
generacies, is explained in more detail in Wuyts et al. 
(2010). We use Hyperz (Bolzonella et al. 2000) to per- 
form SED fitting at a fixed spectroscopic redshift. The 
newest Bruznal & Chariot population synthesis mod- 

http: / / ssc.spitzer.caltech.edu/irac /iracinstrument handbook/ 

4 http:/ /ssc. spitzer. caltech.edu/ mips/ mipsinstrument handbook/ 


els (CB07, kindly made available by the authors - see 
Bruzual & Chariot (2003)) are used, with a Chabrier ini- 
tial mass function (Chabrier 2003) and a Calzetti dust 
extinction law (Calzetti et al. 2000). To reduce the num- 
ber of degrees of freedom and accompanying degenera- 
cies, we fix the metallicity at 0.4 Z 0 for all 4 sources, 
consistent with the abundances indicated by rest-frame 
optical spectroscopy (see §4.1). SED templates in the 
literature most often use constant star formation mod- 
els or exponentially declining models - SFR ~ %~ t l T - 
with a wide range of values for the e-folding time r. The 
current SFR reported by the SED fitting procedure de- 
pends strongly on the assumed star formation history. 
Other stellar population parameters like dust extinction, 
age and even stellar mass can be influenced as well; they 
all work together to produce the observed SED and nu- 
merous degeneracies and trade-offs exist. The range of 
SFHs allowed in the SED fit is therefore not an arbritrary 
choice. Exponentially declining models assume that we 
observe the galaxy at its minimum SFR, which is not nec- 
essarily justified for star- forming galaxies at z~2. Maras- 
ton et al. (2010) explore inverted r-models where the 
star formation rate rises exponentially with time - SFR 
~ e+C r - and find a significant improvement in the re- 
production of stellar population parameters for a sample 
of mock star- forming galaxies at z ~ 2. For this sam- 
ple of lensed galaxies, we choose to limit ourselves to a 
constant star formation history (CSF) as a reasonable 
average over the galaxy’s lifetime and a compromise be- 
tween exponentially declining and rising SFHs. Erb et al, 
(2006c) also find the current SFR to be an adequate rep- 
resentation of the past average SFR for their sample of 
UV-seleeted galaxies at z ~ 2. Additionally, CSF models 
allow a more robust comparison to other SFR indicators 
discussed in §3.4, for which the conversions from luminos- 
ity are generally based on constant star formation stellar 
population models. 

Now that we have fully defined the input SED models, 
we have to consider the allowed range of stellar popula- 
tion parameters, specifically the stellar age. The domi- 
nant emission of O and B stars in the observed SED of a 
star-forming galaxy can cause the inferred stellar pop- 
ulation parameters to be luminosity- weighted towards 
very young models, with unrealistically high SFRs and 
low stellar masses. An often used measure to avoid this 
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luminosity-weighted bias restricts the age of the stellar 
population to be larger than the dynamical timescale of 
the galaxy (Wuyts et ah 2007; Maraston et al. 2010: 
Reddy et ah 2010). From velocity dispersion and size 
measurements of z ~ 2 LBGs, this is inferred to be 
~ 70Myr (Erb et ah 2006c). We illustrate this bias 
with histograms of the best-fit ages and star formation 
rates reported for 1000 mock realizations of the observed 
SEDs consistent with the photometric uncertainties (Fig- 
ure 4). There is a clear bimodality of very young models 
(< 70 Myr) and moderately older models (70-200 Myr), 
which is correlated with a bimodality in the SFR, where 
the young models require a much higher current SFR 
to build up a similar stellar mass. Restricting the age 
to be larger than the dynamical timescale avoids this 
class of unphysically young stellar populations with ex- 
treme SFRs. For cB58, only very young models (ages < 
20 Myr) are found for the mock SEDs, due to the lack 
of a significant Balmer break between the J and H ~ band 
photometry from Ellingson et al. (1996). This was also 
found by Siana et al. (2008), who report a best fit age 
of 9.3i3;[ Myr. However, when the age is restricted to 

> 70 Myr, more than 90 % of the mock SEDs return 
a best- fit model with y 2 < 5.0, which is not unreason- 
able. We will also show in §4.4 that the extreme SFRs 
which accompany very young best-fit models for cB58 
(300-600 Mq yr” 1 as can be seen from Figure 4) com- 
pletely disagree with the other SFR indicators. 

A further argument in favor of the age restriction 
comes from the metallicity. When the metallicity is al- 
lowed to vary between 0.2 Z 0 , 0.4 Z 0 and Z 0 , the SED 
fitting procedure favors a metallicity of 0.2 Z 0 for all 
4 galaxies, lower than indicated by observations of the 
rest-frame optical spectra. When the age is forced to be 

> 70 Myr, a metallicity of 0.4 Z 0 is favored, consistent 
with the results from rest- frame optical spectroscopy. 

Restricting the age of the stellar population limits the 
sensitivity of the SED fit to the most recently formed 
stars (on < 70 Myr timescales) and returns more phys- 
ically meaningful stellar population parameters. A re- 
lated issue is the presence of an old underlying stel- 
lar population (1-2 Gyr timescales), which can contain 
the majority of a galaxy's stellar mass without emitting 
enough light to be detected in the observed SED at rest- 
frame UV and optical wavelengths. The availability of 
IRAC data at rest-frame near-IR wavelengths greatly in- 
creases the leverage, but some contribution of old stars 
cannot be ruled out from single component SED fitting. 
Full-scale multi-component SED fitting can address the 
presence of older stellar populations, but it introduces 
additional degrees of freedom (the number, relative tim- 
ing and relative strength of the different star formation 
episodes), which the data fail to constrain uniquely. We 
can set an upper limit on the contribution of an old stellar 
population to the total stellar mass of the galaxies with 
the extreme case of a very young model, minimally red- 
dened to match the UV spectral slope, combined with 
a maximally old underlying model (Daddi et al 2004, 
Shapley et al. 2005, Wuyts et ah 2010). Specifically, 
we scale a 10 Myr CSF model to the observed r-band 
magnitude, subtract this model from the observed SED, 
and scale the maximally old model (t = 3 Gyr for RC- 
SGA0327 and t - 2 Gyr for SGAS1527, SGAS1226 and 



50 100 150 200 250 

Age [Myr] 



current SFR [Mg/yr] 

Fig. 4. Histograms of the best-fit age (top) and SFR (bottom) 
for 1000 mock SEDs consistent with the photometric uncertain- 
ties. A clear bimodality can be seen. The shaded regions show the 
best- fit SED models with ages < 70 Myr. which correspond exclu- 
sively to the high SFR solutions. This part of parameter space is 
eliminated when the age is restricted to > 70 Myr. 


cB58) to match the residual magnitude at 3.6 ^m. In 
this scenario, the stellar mass is largely dominated by the 
older stellar population, but we find it to be at maximum 
a factor of 1.7 higher than the single-component mass for 
all sources (we have corrected the mass of the older stellar 
population for the mass returned to the ISM by super- 
novae, this is negligible for ages ~1 0-1 00 Myr (Renzini & 
Ciotti 1993)). This places a conservative constraint on 
the maximum contribution of old stellar populations to 
the total stellar mass in these galaxies. 

3.3. Reddening diagnostics 

The presence of dust in galaxies causes reddening. 
Since the stars and the ionized gas are not necessarily 
found in the same locations within a galaxy, the stel- 
lar light and the nebular emission lines may not experi- 
ence the same reddening. In local star- forming galaxies, 
Calzetti et ah (1994) found the reddening of the nebu- 
lar lines to be a factor of ~ 2 higher than the reddening 
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of the stellar continuum. At z ~ 2, the same effect has 
been seen for some star- forming galaxies, both lensed and 
non-lensed (Forster Sehreiber et al. 2009; Finkelstein et 
al. 2009: Bian et al. 2010), but other galaxies show no 
reddening difference between the stars and ionized gas 
(Erb et al. 2006c; Hainline et al. 2009). 

The SED fitting procedure reports E(B — V) S , the red- 
dening of the stellar light. We add the subscript to differ- 
entiate with the reddening experienced by the galaxy’s 
ionized gas, which we will denote as E(B — V) g . This 
reddening can be derived from the Balmer decrement, 
the ratio of the H a to H/3 flux. We assume an intrinsic 
Balmer decrement of 2.85, obtained for recombination 
case B, electron density n e < 10 4 cm 3 and temperature 
T ~ 10 4 K (Osterbrock 1989). When no measure of the 
Balmer decrement is available, we use the empirical rela- 
tion between a galaxy’s stellar mass and its gas reddening 
derived at z ~ 0.1 from the Stripe 82 subsample of the 
Sloan Digital Sky Survey (Eq. 9, Gilbank et. al. 2010). 
This G10 relation has an uncertainty of ~ 0.2 dex on 
E(B — V) g at the relevant stellar masses 5 . This em- 
pirical relation has not been verified directly at higher 
redshift, although agreements with 24 pm emission have 
been found at z ~ 1 (Gilbank et al. 2010), but it is clearly 
preferable to assuming an intrinsic Balmer decrement. 

3.4. Star formation rate diagnostics 

A galaxy’s star formation rate can be estimated from 
several different indicators, the most common of which 
are the UV continuum emission, sometimes combined 
with the infrared emission from dust grains heated by ab- 
sorbed UV photons, and the nebular emission lines from 
the ionized gas. Rest-frame UV luminosity is straight- 
forward to measure at z ~ 2, but the conversion to a 
SFR depends on the galaxy’s stellar population and the 
luminosity requires a large correction for dust obscura- 
tion. Dust corrections can be avoided by including the 
galaxy’s infrared emission, which traces the stellar light 
absorbed by dust grains, such that the combined UV and 
IR luminosity captures the total bolometric output from 
young stars. 

Rest-frame optical emission line indicators require a 
smaller dust correction than the UV continuum, but need 
to be extrapolated from the regions of the galaxy sampled 
by the spectroscopic aperture. These aperture correc- 
tions are especially important for lensed galaxies, where 
the arcs are often too extended to be covered by the 
slit and a precise magnification model is required to ex- 
trapolate the observed line emission to the total intrinsic 
emission of the galaxy. Additionally, in the presence of 
strong stellar population gradients, the observed region 
might not be representative of the galaxy as a whole. 

Star formation from UV luminosity 

Star formation rates based on the UV continuum are 
easy to obtain at z ~ 2. Two or more broadband mag- 
nitudes at rest-frame UV wavelengths suffice to measure 
the UV luminosity and to estimate the dust correction 
from the UV spectral slope, as established by Meurer 
et al. (1999) for local starburst galaxies (based on the 

° Gilbank et al. (2010) use the Kroupa IMF, we convert to the 
Chabrier IMF with !og(M*) C hab = log(M*) K rou P a ” 004 


Calzetti extinction law). However, this method involves 
considerable uncertainties. The conversion between UV 
luminosity and SFR depends strongly on the assumed 
stellar population model, mostly its age and star forma- 
tion history, causing an uncertainty of at least 0.3 dex 
in the calibration (Kennicutt 1998). Additionally, esti- 
mates of the dust content based on the UV spectral slope 
suffer from the degeneracy between a galaxy’s dust con- 
tent, its metallicity and the age of its stellar population, 
where more dust, more metals or older stars will similarly 
cause a redder spectral slope. The availability of a spec- 
troscopic measurement of the metallicity, combined with 
multi-wavelength SED fitting with coverage of the age- 
sensitive 4000A break addresses both problems by fully 
characterizing the galaxy’s stellar population, including 
reliable independent estimates of the dust content, metal- 
licity and age. For these reasons, we will only use the 
SFR reported by the SED fit in our final comparison of 
SFR indicators in §4.4. However, since multi- wavelength 
photometry extending redward of the 4000A break is not 
straightforward to obtain and metallicity measurements 
become increasingly uncertain for galaxies at z > 1, we 
also quantify how reliably we can estimate the dust ex- 
tinction from the UV spectral slope in §4.2. 

Star formation from infrared emission . 

An independent estimate of the dust extinction can 
be obtained from a galaxy’s infrared luminosity. Dust 
grains will re-radiate the absorbed UV emission at in- 
frared wavelengths, such that a galaxy’s combined UV 
and IR luminosity is a good proxy for the total bolomet- 
ric output from its newly formed stars. One caveat is the 
contribution to the infrared emission from dust heated by 
old stars. We can reasonably assume this contribution to 
be negligible for this sample of galaxies, where the ob- 
served SEDs are best fitted by young stellar populations 
and we found strict limits on the presence of older stars 
(see §3.2). We use the SFR conversion from Bell et al. 
(2005) 6 , which is derived for a ^100 Myr old stellar pop- 
ulation with constant star formation. 

M " yr- 1 = 224 X 10_44 ^ /h + l-9iieoo}erg s -1 (1) 

Lift Is the total infrared luminsosity (5-1000 pm) and 
Z /1600 = ^A,i 6 oo* The factor 1.9 accounts for the spec- 
tral slope of the stellar population, such that 1.9Li6oo 
includes all the UV emission. 

Em is preferably evaluated from sufficient infrared and 
submillimeter data to capture the bulk of the bolometric 
energy output of the infrared SED. Unfortunately, data 
at these longer wavelengths is hard to obtain, especially 
for sources at higher redshift, due to confusion noise and 
sensitivity limitations. For this sample, 24 pm data exist 
for RCSGA0327, cB58 and SGAS1527. cB58 has ad- 
ditionally been detected at 70 pm (Siana et al. 2008), 
850 pm (van der Werf et al. 2001) and 1200 pm (Baker 
et al. 2001). Much work has been done to calibrate Ljr 
from observed 24 pm photometry, which traces the mid- 
IR region (6-12 pm) for galaxies at z ~ 2, This wave- 
length range is dominated by polycyclic aromatic hydro- 
carbon (PAH) features, which are stochastically heated 

6 Bell et al. (2005) use the Kroupa IMF, we convert to the 
Chabrier IMF with SFRchab = 0.88 SFR-Kroupa- 
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by single UV photons. Recent studies have shown that 
the physical conditions in star-forming galaxies at z ~ 2 
are similar to those in local star-forming galaxies of sig- 
nificantly lower infrared luminosity (Papovich et ah 2007; 
Rigby et ah 2008; Rujopakarn et ah 2011a), such that 
when local infrared SEDs are used to extrapolate the 
observed 24 pm emission of z ~ 2 galaxies, Ljr will be 
overpredicted by factors of 5-10. This is commonly re- 
ferred to as the mid-IR excess and has been confirmed by 
recent Herschel results (Elbaz et ah 2010; Nor don et ah 
2010). We avoid this bias by using a new prescription for 
L/i?, based on the local IR SED templates from Rieke et 
al. (2009), adapted to z ~ 2 by accounting for the chang- 
ing physical conditions of higher redshift star-forming re- 
gions (Rujopakarn et ah 2011b). This prescription has a 
systematic uncertainty of 0.1 dex. 

Star formation from nebular emission lines . 

The luminosity of the Ha recombination line is directly 
coupled to the incident number of Lyman continuum 
photons produced by young stars, and hence is propor- 
tional to the SFR. Star formation rates derived from dust 
corrected Ha emission are generally seen as robust and 
often used as a comparison and/or calibrator for other 
indicators (e.g. Gilbank et ah 2010, Reddy et al. 2010). 
When the Ha line is redshifted out of the optical win- 
dow, the [O II] A3727 emission line is commonly used as 
an indicator of star formation. The luminosity of forbid- 
den lines is not directly coupled to the number of ionizing 
photons, since the excitation depends on the abundance 
and ionization state of the gas. Gilbank et al. (2010) 
derived a mass-dependent correction for these effects by 
calibrating [O II] SFRs against dust-corrected Ha SFR s 
locally. 

We use the conversions from L Hcx and L\ Q //] to SFR 
from Kennicutt (1998) 7 . The reddening of the nebular 
gas E(B — V) g as derived in §3.3 is used to correct for 
extinction. For [O II], we also present values of the cor- 
rected [O II] SFR following Gilbank et al. (2010) 8 . The 
size of this correction is relatively small over the mass 
range of the galaxies, 

4. RESULTS 

4.1. Physical conditions from i rest- frame optical 
spectroscopy 

Analysis of the nebular emission lines from rest-frame 
optical spectroscopy may allow constraints to be placed 
on extinction, redshift, velocity width, electron density, 
ionization parameter and metallicity of galaxies. The 
results for SGAS1527 and SGAS1226 from the spectra 
presented in §2.2 are summarized below. 

® Extinction 

We detect two Balmer lines for SGAS1527, Ha and 
H 8, in two different filters. Unfortunately, because 
the N7 observation lacks a contemporary telluric 
standard, we do not trust the relative fluxing be- 
tween Ha and H 8 to measure the reddening. Hy 
is not covered, and US is lost to a skyline. For 

* Kennicutt (1998) uses the Saipeter IMF, we convert to the 
Chabrier IMF with SFR Sa ip - 1.7 SFR Chab . ' 

8 Gilbank et ah (2010) use the Kroupa IMF, we convert to the 
Chabrier IMF with SFRchah = 0.88 SFRio-oupa- 


SGAS1226 we detect only one Balmer line, H.0. Hy 
and H<5 were covered by the Nirspec~5b filter, but 
the lines were not detected. Ha is not accessible 
from the ground. Thus, we are not able to reli- 
ably measure an extinction from the Balmer lines 
for either galaxy. 

• Redshift 

We fit the nebular redshift of SGAS1527 using 
the following emission lines: [Ne III] A3869, H/T 
[O III] A4959.5007 and Ha. The derived red- 
shift is z = 2.76195 ± 0.0002. The nebular red- 
shift of SG AS 1226 is derived from the emission 
lines of 113 and [O III] A4959,5007 and yields 
2 = 2.9257 ± 0.0004. Both redshifts are slightly 
higher than those reported in Koester et al (2010), 
which were based on absorption features. The pres- 
ence of galaxy outflows of ~ 200 km/s rest-frame 
can explain these offsets. The emission line red- 
shifts measured here should serve as the systemic 
redshifts, since one expects the H II regions to be 
at rest with respect to the stars. 

• Velocity Width 

We measure velocity widths of bright lines, and 
compare them to the instrumental resolution as 
measured from Argon lamp exposures (see Rigby 
et al, (2011) for more details). Fits to the 
[0 III] A4959,5007 and Ha lines of SGAS1527 result 
in a — 47 ± 4 km s'" 1 . For SGAS1226, a = 100 ± 
20 km s _1 is derived from the [O III] A4959.5007 
lines. 110 and the [O II] A3727 doublet do not 
inform this analysis, since H/3 is too weak for a 
useful measurement of its linewidth, and since the 
[O II] A3727 doublet is not sufficiently resolved for 
a non- degenerate fit. These values are consistent 
with the broad range of a ~ 40 — 200 km s” 1 found 
at z ~ 2 — 3 (Erb et al. 2006b). 

• Electron density 

Rigby et al. (2011) derived an electron density 
from a two-component fit to the flux ratio of the 
[O III A3727 doublet of RCSGA0327. We have at- 
tempted this same procedure for SG AS 1527 and 
SG AS 1226, but the lower signal-to- noise of the 
spectra leads to ambiguous fitting results. In 
SGAS1527, the line profile appears to have two 
components, but fitting it as such returns unphysi- 
cally narrow linewidths - half the instrumental res- 
olution at this wavelength. The source of the prob- 
lem may be the sky line that falls on the bluer 
line of the doublet. As a result, we cannot si- 
multaneously determine both the doublet ratio and 
the velocity dispersion a in the [O II] A3727 dou- 
blet. Instead, we vary a within the range permit- 
ted by the fit to the [O III] A4959.5007 lines, fit- 
ting a doublet ratio for each value. This results 
in a flux ratio £(3727/3729) - 1.03 ± 0.2. Using 
the IRAF task stsdas. analysis. nebular. temden, our 
measurement corresponds to an electron density 
n e = dOOlllo cm ~ 3 at T e = 10 4 K. Following a 
similar procedure for SGAS1226, we find a flux ra- 
tio of f(3727/3729) = 0.59 ±0.1. This corresponds 
to the low-density limit for T e — 10 4 K. 
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Rigby et al. (2011) found a tight constraint of n e = 
235±1| cm” 3 (at T e - 10 4 K) for RCSGA0327. 
They summarize constraints from the literature for 
the Clone, the Cosmic Horseshoe and J 0900+2234 
(see Rigby et ah (2011) for references) as a wide 
range of n e = 600 — 5000 cm" 3 , with large uncer- 
tainties for individual measurements. These con- 
straints are remarkably higher than the low elec- 
tron densities found for RCSGA0327, SGAS1527 
and SGAS1226, which lie closer to the typical den- 
sities of local H II regions, n e ~ 100 cm" 3 (e.g. 
Zaritsky et al. 1994). 

The inverse relation between electron density and 
size of the H II regions found locally (Kim & 
Koo 2001) predicts that local H II regions with 
n e ~ 100 cm” 3 should have diameters of ~ 8 pc. 
If this relation holds at z ~ 2, n e ~ 300 cm” 3 as 
measured on average for RCSGA0327, SGAS1527 
and SGAS1226 corresponds to star-forming regions 
of ~ 3 pc, while n e ~ 1000 cm” 3 as an average for 
the Clone, the Cosmic Horseshoe and J0900+2234 
implies H II regions as small as 1 pc. However, the 
uncertainties on the electron density measurements 
of these latter sources are large; rest-frame optical 
spectra of similar quality to what we present in this 
paper and have published for RCSGA0327 (Rigby 
et al. 2011) is needed for a larger sample of indi- 
vidual star- forming galaxies at z ~ 2. 

• Ionization parameter 

Figure 1 of Kewley & Dopita (2002) illustrates the 
use of the flux ratio of [O III] A5007 to [O II] A3727 
as a diagnostic of the ionization parameter. Assum- 
ing an oxygen abundance of 20-40% solar (on the 
Asplund system), Equation 12 of Kewley & Dopita 
(2002) 9 yields an ionization parameter of log U — 
—2.9 to -2.7 for SGAS1527 and log 17 - -2.8 
to —2.9 for SGAS1226. This assumes no extinc- 
tion; an extinction correction lowers the inferred 
log U. These results increase the number of leased 
galaxies for which this diagnostic is measured from 
four to six. As noted in Rigby et al. (2011), 
they all have very similar ionization paremeters of 
logU = —2.9 to —2.7. This is roughly twice as 
high as —3.2 < log U < —2.9 measured for lower 
luminosity local galaxies (Moustakas et al. 2006. 
2010 ). 

• Oxygen abundance 

Unfortunately, the u gold standard” diagnostic of 
oxygen abundance, 10 III] A4363, is unavailable 
in either galaxy, due to limited wavelength cov- 
erage in SG AS 1527 and skyline contamination in 
SG AS 1226. The ratio of [N II! A6583 to Ha, known 
as the N2 index, is a commonly- used, reddening- 
free abundance indicator. Unfortunately, the reel- 
shift of SGAS1226 is too high for these lines to be 
visible in the near-IR windows. The lines are cov- 
ered in SG AS 1527, though [N II] is not detected, 
limiting us to an upper limit on the oxygen abun- 
dance, We simultaneously fit Ha and the two [N II] 

9 Kewley &: Dopita (2002) use the abundances from Anders Sz 
Grevesse (1989); we convert to Asplund et al. (2009). 


lines, setting a common linewidth for all three lines 
that is allowed to vary within the range measured 
above. The 3 cr upper limit is log ([N II] A6583/Ha) 

< —0.65, which corresponds to an oxygen abun- 
dance of 12 + lo g(0/H) < 8.5 and a metallicity of 

< 70% of solar on the Asplund et al (2009) scale. 
The la upper limit on the oxygen abundance is 
12 + log(0/i7) < 8.3, which translates to < 40% of 
solar. 

The Ne302 index, the ratio of [Ne III] A3869 to 
[O II] A3727 is another accessible, reddening- free 
abundance indicator. Unfort unat ley the [Ne III] 
line is lost to a skyline in the spectrum of 
SGAS1226. For SGAS1527, the measured ratio 
of log ([Ne III] A3869/[0 II] A3727)- — 0.42lg;^f 
yields an abundance of 12 + log(0/i7) = 7.6 dt 0.2 
via the calibration of Shi et al. (2007). This cor- 
responds to a metallicity of 5-12% of solar on the 
Asplund et al. (2009) scale. Unfortunately, as Fig- 
ure 1 of Shi et al. (2007) makes clear, the scatter in 
this calibration is large, 0.3 dex overall and worse 
at low ratios of [Ne III]/[0 II] . This diagnostic does 
therefore not prove useful in this regime. 

4,2. Stellar population parameters 

The stellar populations of the lensed galaxies are con- 
strained by the rest-frame UV to near-IR spectral energy 
distributions. Figure 5 shows the best-fit SEDs for the 
arc and counter-image of RCSGA0327, eB58, SGAS1527 
and SGAS1226 at rest-frame wavelengths. We have in- 
vestigated whether the small discrepancy between the 
photometry and the best-fit SED model that can be seen 
for RCSGA0327, cB58 and SGAS1226 around the 4000i 
break can be significantly reduced by allowing exponen- 
tially declining SFHs. We find that this is not the case 
for cB5S and SG AS 1226. For these galaxies the discrep- 
ancy is due to a tension between our minimum age con- 
straint of 70Myr (well motivated in §3.2) and the lack 
of a significant 4000A break in the observed photometry. 
For RGSGA0327, SED models with an e-folding time of 
r = 10 — 50Myr provide improved fits to the photometry, 
shown by the dashed green lines in Figure 5. The current 
SFR derived from the SED fit depends greatly on the as- 
sumed star formation history. For both the arc and the 
counter-image of RCSGA0327, exponentially declining 
SFHs with r = 10 — 50 My r produce SFRs < 20 M 0 yr” 1 , 
much smaller than indicated by the other SFR measure- 
ments (see §4.4). Recently, Wuyts et al. (2011) have sim- 
ilarly reported the need for a minimum e-folding time of 
log(r) = 9.5 (for CB07models) to avoid SFRs that are 
underestimated compared to measurements derived from 
the bolometric luminosity for galaxies out to z ^ 3. For 
the relatively young ages of the lensed galaxies in this 
sample, such long e- folding times are equivalent to the 
assumption of constant star formation. 

Table 3 reports the stellar population parameters de- 
rived from the SED fit. The confidence intervals are 
estimated from 1000 mock realizations of the observed 
SEDs which are consistent with the photometric uncer- 
tainties, and by varying the age restriction between 50 
and 100 Myr. The age t and reddening factor E(B — V) 
depend solely on the galaxy’s colors; the stellar mass 
and star formation rate have an additional dependence 
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Fig. 5. Best fit stellar energy distributions for the arc (top) 
and counter-image (bottom) of RCSGA0327, cB58, SGAS1527 and 
SGAS1226 derived from age-restricted CSF models. The dashed 
green lines show the improved SED fits for both the arc and 
counter-image of RCSGA0327 when using an exponentially declin- 
ing SFH with e-folding time r=10-50 Myr. Photometry datapoints 
are overplotted with 1 o errorbars. 


on the lensing magnification. Based on WFC3/HST 
images, Sharon et al. (2011) (in preparation) report a 
magnification factor of 3.0±0.2 for the counter- image of 
RCSGA0327 and an average magnification of 25.1 ±3.2 
across the giant arc. Koester et al. (2010) report mag- 
nification limits of < 15 for SGAS1527 and > 40 for 
SGAS1226. cB58 is magnified by a factor of 30, with an 
estimated uncertainty of 30% (Seitz et al. 1998). The 
magnification uncertainties are included in the reported 
uncertainties for the stellar mass and SFR. 

The stellar populations of all four lensed galaxies are 
remarkably similar and can be characterized on average 
with an age ~ 110 Myr, reddening E(B — V) s ~ 0.17, 
SFR ~ 65 M@ yr” 1 and stellar mass ~ 5.5 x 10 9 Mq. 
Compared to representative samples of U V -selected star- 
forming galaxies at z ~ 2 — 3 (Shapley et al. 2001; Erb et 
al. 2006b; Mannucci et al. 2009: Reddy et al. 2010), these 
lensed galaxies represent young, star bursting systems 
with low stellar masses, suggesting they have only re- 
cently begun building their stellar population. The same 
conclusion is reached from a comparison of the UV lumi- 
nosity and stellar mass to the characteristic luminosity 
and stellar mass of comparable galaxies at z ~ 2 — 3. Us- 
ing the characteristic luminosity from Oesch et al. (2010) 
at 0 — 1.5 — 2.0 for a comparison to RCSGA0327 and the 
results from Reddy & Steidei (2009) at z — 2.7 — 3.4 for 
a comparison to cB5S, SGAS1527 and SGAS1226, we 
find an average brightness of ~ 2.6 L*. The stellar mass 
of the sample lies at ~ 0.08 M* tar when comparing to 
the characteristic stellar mass found by Marehesini et al. 
(2009) at z ^ 2 — 3. 

We now revisit the issue raised in §3.4 of how reliably 
one can estimate the dust extinction E(B — V) s from 
the UV spectral slope alone. This is relevant for galaxies 
at z > 1, where the photometric and spectroscopic data 
required to fully characterize the galaxy’s stellar popula- 
tion through SED fitting becomes hard to obtain. We use 


the UV color of the galaxies, derived from the filters that 
best match the rest-frame 1600A-2300A window, to mea- 
sure the UV spectral slope 3 (f\ ~ A 5 ). The reddening 
E(B~ V) s derived from 0 and the relation established by 
Meurer et al. (1999) for local starburst galaxies is consis- 
tently a factor of 1.25-2.0 lower than the value reported 
by the SED fit. This underestimates the dust-corrected 
SFR of the galaxies by factors of 2-3.5. 

We investigate the origin of this discrepancy by es- 
tablishing the relation between 3 and E(B — V) s for 
representative CB07 models. For a range of values 
E(B — V) $ — [0,0.3] we redden each SED model with 
the Calzetti dust extinction law and measure the UV 
spectral slope based on the same filters available to us 
for each of the four galaxies in the sample (specifically 
B and r for RCSGA0327, V and I for cB58, r and z 
for SGAS1527 and r and i for SGAS1226). Since 0 de- 
pends somehat on the exact wavelength window used, 
this results in four slightly different relations between 3 
and E(B ~~ V) s for each SED model. Figure 6 shows the 
Meurer relation in red, offset from the E(B — values 
reported from the SED fit (black datapoints). The rela- 
tions between 0 and E(B — V) s for a solar metallicity, 
500 Myr old CSF model and a 100 Myr old CSF model 
with a metallicity of and 0.4 Z@ are shown in cyan, 
blue and green respectively. The 100 Myr, 0.4 Zq, CSF 
model is representative of the average stellar population 
of the galaxies, and can be seen to agree with the best- 
fit values for E(B — V) 8 from the SED fit. The lower 
metallicity accounts for much of the discrepancy with 
the Meurer relation: a 0.4 Zq SED model is intrinsically 
less red than a solar metallicity SED model and thus re- 
quires a larger dust correction to match the same UV 
spectral slope. In general, this highlights the danger of 
blindly adopting the Meurer relation to dust-correct the 
UV continuum emission of high redshift galaxies, with- 
out additional information on the age and metallicity of 
the stellar populations. 



“ 2.4 “ 2.2 - 2.0 - 1.8 - 1.6 -* 1.4 - 1.2 - 1.0 

P 

FIG. 6. The relation between the UV spectral slope / 3 and the 
reddening E(B — V) s - The local relation established by Meurer et 
al, (1999) is shown in red. The black datapoints show the galaxies' 
best- fit E(B — V)s from the SED fit. The relations we derived for 
the SED models are shown in cyan, blue ana green. The 100 Myr, 
0.4 Z©, CSF model is representative of the average stellar popula- 
tion of the galaxies and agrees with the results from the SED fit. 
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Table 3. Stellar population parameters. 


source 

Age 

(Myr) 

E(B - V) s 

log(M»/M®) 

SFR 

(M 0 yr -1 ) 

RCSGA0327-A 

128±13 

0.15±0.035 

9.83±0.06 

71±12 

RCSGA0327-C 

143 ±42 

0.10±0.04 

9.76±0.05 

52+15 

cB58 

72±18 

0,25±0.055 

9.69±0.14 

82 ±2 7 

SGAS1527 

126±31 

0.15+0.05 

> 9.82 ± 0.04 

> 71 ± 16 

SGAS1226 

86±41 

G.20±0.05 

< 9.54 ±0.07 

< 48 ± 15 


Table 4. Extinction estimates. 


source 

E(B - V) s 

E(B - V r )^ aZmer 

E(B - V)g- 

RCSGA0327-A 

0.15±0.035 

0.20±0.19 

0.31±0.18 

RCSGA0327-C 

0.10±0.04 


0.29T0.17 

cB58 

0.25±0.055 

O.lliO.lO 

0.2Q±0.12 

SGAS1527 

0.15±0.05 


0.24±0.14 

SGAS1226 

0.20±0.05 


0.15±0.09 


Column 1: Reddening of the stellar light from the SED fit 
Column 2: Reddening of the nebular gas from the Balmer decre- 
ment 

Column 3: Reddening of the nebular gas based on the stellar 
mass as derived by Gilbank et al. (2010) 


4.3. Reddening 

We report the reddening of the stellar light, E(B — V) S 
and the reddening of the ionized gas, E(B — V) g in Ta- 
ble 4. These measures of reddening are important to 
reliably correct the SFR indicators for dust extinction. 
E(B — V) s is taken from the SED fit. E(B ~~ V) g is es- 
timated from the stellar mass following the local G10 
relation (Gilbank et al. 2010, see §3.3). For the giant arc 
of RCSGA0327 and for cB58, a measurement of the gas 
reddening can also be made from the Balmer decrement 
(Rigby et al. 2011; Teplitz et al. 2000) and the indepen- 
dent estimates of E(B — V) g appear consistent within 
the la uncertainties. This suggests that the local G10 
relation remains valid at z ~ 2 and there is mild or no 
evolution in the dust content of galaxies as a function 
of mass between these epochs. However, we note that 
this is only based on two measurements and more work 
is needed to verify the validity of this local prescription 
at higher redshift. For RCSGA0327 and cB58, we use 
the weighted average of both measures of E(B ~~ V) g 
to correct the nebular emission line SFR indicators for 
extinction. 

Within the la uncertainties, these reddening measure- 
ments cannot conclude whether or not the nebular emis- 
sion lines experience larger reddening than the stellar 
light in this sample of star- forming galaxies. 

4.4. Star Formation Rates 

We have outlined 4 different SFR indicators in §3.4 
based on the SED fit, the total bolometric luminosity 
and the Ha and [O II] A3727 emission lines. The SFR 
result from the SED fit is reported in §4.2, after being 
corrected for the lensing magnification (25.1 ± 3.2 for 
RCSGA0327-A, 3.0 ± 0.2 for RCSGA0327-C, 30 ± 9 for 
cB58, < 15 for SGAS1527 and > 40 for SGAS1226). We 
summarize the de-lensed UV and IR luminosities of the 


galaxies (corrected using the same magnification factors) 
in Table 5 and use Equation 1 to derive the SFR from 
the bolometric luminosity. 

The emission line luminosities of SGAS1527 and 
SGAS1226 are derived from the observed fluxes reported 
in Table 2. As can be seen in Figure 1, the N I RS PEC slit 
covers the whole source for both galaxies, so we can use 
the total magnifications of < 15 for SGAS1527 and > 40 
for SGAS1226 to calculate the de-lensed luminosities, 
without additional aperture corrections. Since no mea- 
surement of Ha exists for either source (we don’t trust 
the fluxing of SGAS1527 due to a time delay with the ob- 
servation of the telluric star and for SGAS1226 Ha falls 
outside the observable window due to its high redshift), 
we use 110 and the extinction E(B — V) g derived from 
the stellar mass following the G10 relation (see § 4.3) to 
estimate L# a . For cB58, the Ha and [O II] A3727 lumi- 
nosities published by Teplitz et al. (2000) are used and 
corrected for a total magnification of 30 ± 9. For RC~ 
SGA0327, we use the emission line measurements from 
Rigby et al. (2011). Sharon et al. (2011, in preparation) 
calculate an average magnification of 42.2 ± 5.5 for the 
area of the arc covered by the NIRSPEC slit. Rigby et 
al. (2011) recommend using the observed Up flux and the 
Balmer decrement to estimate L# a . For all galaxies, de- 
lensed Ha and [O II] A3727 luminosities are summarized 
in Table 5 and the resulting SFRs are dust-corrected us- 
ing the Calzetti extinction taw and E(B — V) g from § 4.3. 

The different SFR indicators are compiled in Table 6 
and Figure 7 for easy comparison. For SG AS 1527 and 
SG AS 1226, all SFR measurements are upper/lower lim- 
its due to the poorly constrained magnification; however 
on a relative scale, these limits can still be used for a 
comparison of the different indicators. 

The complex lensing signature of RCSGA0327 does not 
allow a direct comparison of the SFRs between the gi- 
ant arc and the counter-image. The lens model tells us 
that the counter-image is a more or less uniformly mag- 
nified image of the source plane galaxy. The giant arc 
on the other hand consists of three merged images of the 
source galaxy and depending on the location with re- 
spect to the lensing caustic, different regions within the 
source plane galaxy receive very different magnifications 
(see Sharon et al. 2011, in preparation, for more details). 
The de-lensed luminosities and derived SFRs of the giant 
arc are therefore not representative of the source galaxy, 
and should not be compared directly to the results from 
the counter-image. Additionally, Sharon et al. 2011 (in 
preparation) find that the area of the arc observed spec- 
troscopically translates into a very small region of the 
galaxy in the source plane, which consists of a few highly 
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Table 5, De-lensed luminosities. 


source 

L 1600 

10 44 erg s~ l 

Lir 

10 44 erg s“ x 

Lh q 

10 42 erg s~ x 

L[o in 
10 42 erg s '" 1 

RCSGA0327-A 

4.3=fc:0.6 

8.4±3.9 

5.4±1.4 a 

3.3±0.4 a 

RCSGA0327-C 

5.1±0.5 

2.4±1.3 



cB58 

2.1±0.7 

4.9±2.2 

2.6±0.8 b 

2.6±0.8 b 

SGAS1527 

> 5.3 ±0.4 

~3.7±1.7 C 

> 14 ± 8 d 

> 4.6 ±0.9 

SGAS1226 

< 2.0 ±0.1 


< 1.1 ±0.4 d 

< 1.7 ±0.3 


a Rigby et at. (2011) - corrected for an average magnification of 42.2 ± 5.5 as calculated by Sharon et al. (2011, in preparation). These 
luminosities are derived from the area of the giant arc targeted for spectroscopy, which translates into a small region of the galaxy in the 
source plane (Sharon et al. 2011, in preparation), and are not valid for the whole galaxy. 
b Teplitz et al. (2000) 

C SGAS1527 was not detected at 24 /im, the reported luminosity represents a 3 a upper limit, corrected for a lower limit on the magnification 
of < 15 

d Ha is derived from H/3 and the Balmer decrement as derived from the stellar mass through the G10 relation (Gilbank et al. 2010) 


Table 6. Star Formation Rates. 


source 

SFRs ED 

SFRuv+ir 

SFRffa 

SFR [0 w 

SFRseD-SMC 

RCSGA0327-A 

71±12 

37±9 

44±17 a 

48±21 a 

37±6 

RCSGA0327-C 

52±15 

27±4 



39±7 

cB58 

82±27 

20±6 

16±6 

29±13 

38±14 

SGAS1527 

> 71 ± 16 

~ 31 ± 4 b 

> 112 ±73 

> 66 ± 29 

39±4 

SGAS1226 

< 48 ± 15 


< 7 ± 3 

< 20 ±8 

16±2 


a These SFRs are derived from the area of the giant arc targeted for spectroscopy, which translates into a small region of the galaxy in 
the source plane (Sharon et al. 2011, in preparation), and are not valid for the whole galaxy. 
b SGASl527 was not detected at 24 /xmu the reported luminosity represents a 3cr upper limit, corrected for a lower limit on the magnification 
of < 15 


magnified individual star-forming complexes. The SFRs 
measured for this area from the Ha and [O II] A3727 
emission lines are not valid for the whole galaxy and we 
do not include them in our overall comparison of SFR in- 
dicators. A future paper will constrain the Ha luminosity 
of the source galaxy from rest-frame optical spectroscopy 
of the counter- image. 
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Fig. 7. Comparison of SFR, indicators. The gray stars corre- 
spond to the mass-dependent G10 correction for SFR\q jji from 
Gilbank et al. (2010). The SFR, derived from the SED fit is overall 
larger than the other indicators, which suggests that the Calzetti 
extinction law overpredicts the dust correction for the UV luminos- 
ity. The results from an SED fit with the steeper SMC extinction 
law are shown as gray circles. 


Figure 7 shows a general trend where the SFR derived 
from the SED fit is overall larger than the other indica- 
tors (except for SGAS1527). This suggests that the SED 
fit, and thus the dust-corrected UV continuum emission, 
overpredicts the star formation rate. A similar result was 
found for young LBGs (age < 100 Myr) z ~ 2 (Reddy 
et al. 2010), and for cB58 and the Cosmic Eye (Siana et 
al. 2009). These authors have suggested that the Calzetti 
extinction law is not a good representation of the dust 
geometry in young LBGs at z ~ 2 and will significantly 
overpredict the dust extinction. 

The influence of the shape of the extinction law on the 
dust extinction is commonly illustrated with the Meurer- 
plot (Meurer et al. 1999), plotting the UV spectral slope 
8 (fx ~ X p ) versus the ratio of IR to UV luminosity, 
which serves as a parametrization of the dust extinc- 
tion 11 . On this plot, a steeper extinction law results in 
comparatively more absorption of photons at blue wave- 
lenghts and thus a redder UV spectral slope 8 for the 
same extinction. For local star bursts, plotted as little 
black crosses in Figure 8, the Calzetti extinction law is 
the best fit (Meurer et al. 1999). Reddy et al. (2010) 
have shown that at 2 ~ 2, the Calzetti law remains valid 
for star- forming galaxies with ages > 100 Myr. However, 
the younger LBGs in their sample, as well as cB58 and 
the Cosmic Eye, all tend to lie below the Calzetti curve, 
showing less extinction than would be derived from the 
UV spectral slope. For these sources, the steeper extinc- 

11 Meurer et al. (1999) define this relation for Lfir , integrated 
from 40 to 120 where Lir ~ 1.75 x Lfir based on the cali- 
bration by Calzetti et al. (2000). 
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tion law derived for the Small Magellanic Cloud (SMC, 
Prevot et al. 1984) seems to be a better fit. We can add 
datapoints for the arc and counter-image of RCSGAG327, 
for SGAS1527 and for our re-analysis of cB58 to this plot. 
Due to the relatively blue spectral slopes of RCSGA0327 
and SG AS 1527, these galaxies have limited leverage to 
distinguish between the different extinction laws. 



FlG. 8 . Dust attenuation, parameterized as the ratio between 
IR and UV luminosity, versus rest-frame UV spectral slope /?. The 
measurements of local starburst galaxies are plotted as little black 
crosses. The best-fit- relation corresponds to the Calzetti dust ex- 
tinction law (Meurer et al 1999). The steeper SMC extinction law 
is also shown. The red datapoints for cB58 and the Cosmic Eye are 
taken from Siana et al. (2008) and Siana et al. (2009) respectively. 
The relatively blue spectral slopes of RCSGA0327 and SG AS 1527 
have limited leverage to distinguish between the different extinc- 
tion laws. 

We have repeated the SED fitting with the steeper 
SMC extinction law. With the exception of cB58, the 
best-fit SED can not be distinguished from the result 
obtained with the Calzetti law in terms of the reduced 
X 2 of the fit. The improved fit for cB58 was also reported 
by Siana et al. (2008). The derived SFRs are lower by 
a factor of 1.3-2 compared to the Calzetti extinction law 
and therefore in closer general agreement to the other 
SFR indicators (gray circles in Figure 7). With the SMC 
law, the stellar mass is consistent with previous results, 
but the best-fit ages are higher by a factor of ~ 2 and 
there is no sub-population of very young extreme star- 
bursts, which makes the age restriction to older than 
70 Myr unnecessary. 

5. SUMMARY 

This paper presents a uniform analysis of rest-frame 
UV to near-IR spectral energy distributions and rest- 
frame optical spectra for a sample of four of the bright- 
est lensed galaxies at z ~ 2. Given the limitations of 
current facilities, such extensive data collection and de- 
tailed analysis of stellar populations and physical condi- 
tions would not be possible for individual galaxies at this 
redshift without the magnification induced by gravita- 
tional lensing. SED fits to the rest-frame UV to near-IR 
photometry result in very uniform stellar population pa- 


rameters across the sample, which we can characterize as 
young, bright, star bursting systems with little dust con- 
tent and low stellar masses. We derive a strict limit on 
the presence of an older stellar population. Restricting 
the age of the stellar population to be larger than the 
dynamical timescale of 70 Myr during the SED fitting 
is necessary to avoid extreme young starburst models 
and to obtain metallicities consistent with those found 
from the rest-frame optical spectroscopy. When the lo- 
cal Meurer relation (Meurer et al. 1999) is used to es- 
timate the dust correction from the UV spectral slope, 
E(B — V) S is underestimated by factors of 1.25-2, leading 
to underestimates of the dust-corrected UV continuum 
SFR by factors of 2-3.5. This is largely due to the lower 
metallicity (0.4 Z 0 ) of the galaxies. This result cautions 
the use of the Meurer relation without further charac- 
terization of a galaxy’s stellar population and motivates 
efforts to measure metallicities for larger samples of both 
lensed and field galaxies at z ~ 2. 

We present SFR estimates based on the SED fit, 
the UV+IR bolometric luminosity and the Ha and 
[O II] A3727 emission lines. The SED fit with the default 
Calzetti extinction law reports a SFR estimate that is 
consistently higher than the other indicators. Similar re- 
sults were found for cB58 and the Cosmic Eye (Siana et 
al. 2009). This suggests that the Calzetti extinction law 
is too flat for young star- forming galaxies at z ~ 2 and 
therefore over predicts the dust extinction. In the diag- 
nostic Meui'er plot, the galaxies’ spectral slopes are un- 
fortunately too blue to differentiate between the Calzetti 
law and the steeper SMC extinction law. SED fitting 
with the SMC law returns SFRs lower by a factor of 1.3- 
2, bringing this method into general agreement with the 
other SFR indicators, within the la uncertainties. These 
results are a first attempt to compare the common SFR 
indicators in detail for individual galaxies at z ~ 2; the 
comparison would clearly benefit from a larger sample of 
galaxies. 
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